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O . Abstract 



Treating the QCD Wilson loop as amplitude for the propagation of the first 

^ . quantized particle we develop the second quantization of the same propagation. 

The operator of the particle position X^ (the endpoint of the "open string") is 

introduced as a limit of the large A^ Hermitean matrix. We then derive the set 

Q^l of equations for the expectation values of the vertex operators ( V{ki) . . . V{kn) )■ 

The remarkable property of these equations is that they can be expanded at small 
momenta (less than the QCD mass scale), and solved for expansion coefficients. 
This provides the relations for multiple commutators of position operator, which 
can be used to construct this operator. We employ the noncommutative probability 

5r ! theory and find the expansion of the operator ^^ in terms of products of creation 

operators a|j. In general, there are some free parameters left in this expansion. In 
two dimensions we fix parameters uniquely from the symplectic invariance. The 
Fock space of our theory is much smaller than that of perturbative QCD, where the 
creation and annihilation operators were labelled by continuous momenta. In our 
case this is a space generated by d = 4 creation operators. The corresponding states 
are given by all sentences made of the four letter words. We discuss the implication 
of this construction for the mass spectra of mesons and glueballs. 
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1 Introduction 

The interpretation of the Wilson loop W{C) as an amplitude of the charge propagation has 
a long history. Already in early Feynman's work the abelian phase factor was interpreted 
this way. Wilson elaborated on that in his famous paper " Confinement of Quarks" , taking 
the first quantized picture. 

The quark propagator in this picture is given by the sum over closed loops C of the 
loop amplitudes J2c^{C)G{C), where G{C) is the amplitude for the free quark. In 
the large A^ limit there are no dynamical quarks, which makes this picture especially 
attractive. There is a tempting analogy with the string theory. In fact, the string theory 
originated as the planar diagram theory of quark confinement and was boosted by the 't 
Hooft's discovery [|| of the large N limit. 

Now, what would be the second quantized picture of this process? This is not just the 
matter of mathematical curiosity. The development of Feynman- Wilson picture for the 
last 20 years led us to the dead end. We derived^] nonlinear loop equation for W{C), 
but this equation turned out to be too hard to solve. 

In two dimensions, where there are only global degrees of freedom, related to topology 
of the loop, the loop equations were solved^]. In four dimensions the best bet was the 
Eguchi-Kawai reduction, which, unfortunately worked only on the lattice. Taking the 
continuum limit in the EK model appears to be as hard as in the original partition 
function of the lattice theory. This forces us to look for the new approaches, such as the 
second quantized picture. 



1.1 Classical vs Quantum Loops 

The second quantization of the usual field theory involves the following steps 

fields =^ operators =^ commutators =^ H'^ = E^^. (1) 

We would like to follow these steps in QCD for the particle position field X. The para- 
metric invariance of the Wilson loop would make this one dimensional topological field 
theory. The Hamiltonian H vanishes in such theory, as the field operators do not depend 
of the (proper) time. 

To understand this striking phenomenon, let us consider the following simple example: 
Wilson loop in constant abelian electromagnetic field F^jy 

WiC) = exp (^-t F^, £ x^ dx,^ . (2) 

The exponential here represents symplectic form which immediately leads to the commu- 
tation relations 

[X,,X,]=i{F-')^,. (3) 

of \d harmonic oscillators. 

This corresponds to the following momentum space path integral^ 

W{C) oc ^ j Dp{.)5\p{so)) exp {i j dx^p^ tr f exp (i j> dp^X^ , (4) 

One may regard p'As) as Lagrange multiplier for the constraint on the coordinate x^{s) ~ 
X^ modulo total translations 6x^{s) = const. The constraint 

p{so) = (5) 

is needed to eliminate the translations oi p{s). The choice of sq is arbitrary, it is a matter 
of convenience. 

The trace is calculable for the oscillator, it yields 

tr Texp U j dp^xA oc exp [--i{F'^)^y j P^idp^j . (6) 

The integration over p{.) correctly reproduces the initial formula for W{C). For the 
readers reference, this functional integral is computed in Appendix A, including the nor- 
malization (irrelevant for our present purposes). 

Note, that the operators X^ do not commute, and, hence, cannot be simultaneously 
diagonalized. So, we cannot integrate out the Lagrange multiplier and claim that x^{s) 
coincides with one of the eigenvalues of X^. 

^ Later we shall slightly modify this formula, to take into account the discontinuities of the local 
momentum p{s). 



The Wilson loop enters the physical observables inside the path integral, such as 

r dT I Dx{.)6'^{x{sq)) exp l-Tdt (h^ + mA j W{C) (7) 

for the scalar particle propagator. With our representation after integrating out x{.) this 
becomes 

j^ dT j Dp{.)fe^Y> i[dt (i p^X^ - Hip))] (8) 

with the (proper time) Hamiltonian 

H = m'+pl. (9) 

For the Dirac particle we get the same formula with the Dirac Hamiltonian 

H = m + i-ff,p^ (10) 

The only difference with the usual Hamiltonian dynamics is the replacement of the particle 
position x{t) by the operatorX in the px term. The implications of this replacements will 
be studied later. 

1.2 Master Field and Eguchi-Kawai Reduction 

In a way, our approach is an implementation of the Witten's master field ideap. He 
suggested to look for the classical x-independent field A^ in coordinate representation 

W{C) = i^tr texp [j^dx.A^ . (11) 

Clearly, the naive implementation, with A^ satisfying the Yang-Mills equations, does not 
work. There are corrections from the fluctuations of the gauge field. These corrections 
are taken into account in the loop equations[0, |^, so it is worth trying to find effective 
master field equation for A^ from the loop equations. 

In some broad sense the master field always exists. It merely represents the covariant 
derivative operator 

V^ = -zP^ + ^^(0), (12) 

where P^ is momentum operator and A^ is the gauge field operator in the Hilbert space. 
As one may readily check, the matrix trace of the ordered exponential of this operator 
reduces to the operator version of the loop exponential times the unit operator in Hilbert 
space 

tr Texp ( / V ^dx^ j = tr Texp ( / Afj,{x)dx^ j , (13) 

A^{x) = exp {-iVuX^) Af,{0) exp (iVuX^) . 



Taking the trace of this relation in Hilbert space of large A^ QCD, and dividing by 
trace of unity, we would get the usual Wilson loop. We could avoid division by an infinite 
factor by taking the vacuum average instead. 

(0| Texp ^ V^dx^^ |0) = -tr f exp ^ A^{x)dx^^ (14) 

This formal argument does not tell us how to compute this master field. The operator 
version of the Yang-Mills equations is not very helpful. It reads 

N„N„^u]]=9!ff^- (15) 

The operator ^^ in turn, is given by its commutator relations with V^. 

Formally, inserting the left side of this equation inside the trace in the master field 
Anzatz ([Tl|) we would reproduce the correct loop equation provided the trace of the open 
loop is proportional to the delta function 6'^{x — y). The normalization would come out 
wrong, though, unless we include the divergent constant S'^{0) = A'^ in effective coupling 
constant 

glff = gl^' (16) 

The Eguchi-Kawai reduction]^ was, in fact, a realization of this scheme. There, the 
large N matrix A^ represented the covariant derivative operator in above sense, except it 
was not a single matrix. Its diagonal components played the role of momenta A** = — ^p* . 
These were classical, or, to be more precise, quenched^] momenta. 

The off-diagonal components were fluctuating with the weight exp {—A.'^Lym\ corre- 
sponding to the Yang-Mills Lagrangean for the constant matrix field. The cutoff A explic- 
itly entered these formulas. In effect, one could not take the local limit in the quenched 
EK model. Apparently something went wrong with implementation of the master field. 



1.3 Noncommutative Probability Theory 



Recently the large A^ methodology was enriched |T0|, pl]| by using so called noncommutative 



probability theory]^. In particular, in the last paper |]TT| the master field was constructed 
by similar method in the QCD2. 

This amounts to the deformation of the commutation relation 

a^,{p)al{p) - q al{p)a^{p) = 6^,J'^{p - p) (17) 

Instead of the usual Bose commutation relations, with g = 1 one should take q = 0, which 
makes it so called Cuntz algebra. This algebra in the context of the planar diagram theory 
in QCD was first discovered long ago by Cvitanovic et.al.§]. 



As noted in[^, this algebra corresponds to Boltzmann statistics. The origin of the 



Boltzmann statistics in the large A^ limit is clear: at large rank each component of the 



matrix can be treated as unique object, neglecting the probability of the coincidences of 
indexes. The indexes could be dropped after this. 

This interesting construction proves the concept, but we cannot be completely satisfied 
with the fact that the operator representation in this paper is gauge dependent. It used 
the fact that in the axial gauge there is no interaction in QCD2, so that the only deviation 
from the free field behaviour of the Wilson loop comes from the noncommutativity of the 
(Gaussian) gauge potential at various points in space. 

This theory is purely kinematical, it does not make use of the large N dynamics. The 
Fock space is only a slight modification of the perturbative Fock space. The color con- 
finement was not taken in consideration. There should exist a much tighter construction, 
built around the nonperturbative Fock space. 

Also, the gauge dependence should be absent in the final formulation of the theory. 
This resembles the early attempts to solve the Quantum Gravity in 2D by means of 
perturbative Liouville theory. Taking the conformal gauge and using the methods of 
conformal field theory, one was forced to work in the space of free 2D particles, with the 
Bose operators a„, ajj. 

However, as we know now, the Fock space of 2D Gravity is much smaller. It is 
described by the matrix models (or topological field theories), with only discrete states. 
If we now employ the same noncommutative probability theory, we could reinterpret the 
matrix models of 2D Gravity in terms of a single pair of operators a, aK 

To conclude this discussion, the ideas of the noncommutative probability are quite 
appealing, but implementation of these ideas in QCD is not yet finished. The correct 
nonperturbative Fock space in confining phase of QCD is yet to be found. 

1.4 Position Master Field 

In the present work, we follow the same general large N philosophy, but interpret the 
master field as the momentum operator P^ = tA^j, of the endpoint of the QCD string. 
With the second quantization, developed below, we rather construct the position operator 
X. After that we shall reconstruct the momentum operator as well. 

The advantage of the position operator is that (up to irrelevant global translational 
mode) we expect its spectrum to be bounded from above in confining phase. 

The commutators [X^,X,^] describe the "area inside the Wilson loop". We could work 
with expectation values with products of such operators. As for the momentum operators, 
the high momentum region corresponds to small distance singularities of perturbative 
QCD. It makes the matrix elements of commutators of the Witten's master field infinite. 

To rephrase it, the Wilson loop in coordinate space is singular. It cannot be expanded 
in powers of the size of the loop. There are the gluon exchange graphs, which lead to 
singularities at intersections and corners. The gluon exchange graphs for W{C) scale as 



powers of |C|^ '^. 

In two dimensions the singularities reduce to existence of various topological sectors, 
with different Taylor expansions. Without self-intersections, there is a pure area law, with 
intersections there are polynomials of areas inside petals times exponentials]^]. 

In four dimensions there are logarithms, corresponding to the asymptotic freedom. It 
is not clear how to reproduce these properties by any constant large A^ matrix A^. 

With the position operator X, as we shall see, the correct coordinate dependence is 
built in our Anzatz (^. In particular, the 6 function in the loop equation comes out 
automatically, without any assumptions about spectrum of X. We escape the quenched 
EK model paradox of extra factor A'' in effective coupling constant. 

2 Loop Equations 

In this section we review and further develop the method of the loop equations in the 
large A^ QCD, first in coordinate space, then in momentum space. 

2.1 Spatial Loops 

We describe these loops by a periodic function of a proper time s 

C : x^{s) , x^{s + 1) = x^{s). (18) 

The following Schwinger-Dyson equation is to be used 

= lDAtrl'dta'Ma',jj^^ (19) 

exp Ud'x*^^ texp d^' dsA„(x(s))x'^(s)j . 

The variation of the first exponential yields the classical part of equation of motion, 

whereas the second (ordered) exponential provides the commutator terms, present only 
for the self-intersecting loops. 

These terms involve g^ § dti § dt25'^ {^iii) — ^{^2)) ^fiu times the product of two ordered 
exponentials, corresponding to the two petals 6*12,6*21 of our loop C. The S function 
ensures that both these petals are closed loops, as it is required by the gauge invariance. 

^ These self-intersections are responsible for all the "quantum corrections" to the classical equations 
of motion. Dropping them as "rare" or "exceptional" configurations would bring us back to the classical 
Yang-Mills theory. 




The classical part can be represented 0, |T4| as the functional Laplace (or Levy) oper 



ator 



ati 



ii+O § 

dt2 



Sx^{ti) Jti-0 Sx^{t2) 



acting on the loop functional 



W[C] = ^^tTfexpi^£dsA,ixis))x'^is: 



(20) 



f211 



This Levy operator is the only parametric invariant second order operator, made of func- 
tional derivatives. It picks up the contact term in the second functional derivative 



5W[C] 



5x^{ti)6x^{t2) 



oc 



:22) 



5(ti -t2)<(ti) ^tif[D^F^,{x{h))] exp (j^ A^{x{s))x'^{s)\ 
+regular terms. 

On the right hand side, the average product of two loop functionals reduces at large N 
to the product of averages 

— trTexpf / dsiAf,{x{si))x'^{si)] (23) 

— trTexp( / ds2A^,{x{s2))x'^{s-2) 

We thus arrive at the conventional loop equation 

LW\C\ = Ngljdhx'^ih) j'" dt2x'^{t2)5'' {x{t^) - x{t2))W[C^2]W[C2il (24) 

where the second integral goes around the loop avoiding the first point ti. Thus, the 
trivial self-intersection point at ^2 = ^i is eliminated (see0 for details). 

Singularities in this equation should be properly understood. The double integral of 
the d dimensional 6 function does not represent an ordinary number in more than two 
dimensions. This is a singular functional in the loop space, the only meaning of which 
could be provided by an integration with some weight function. 
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Otherwise we should mess up with the gauge invariant regularization of the delta func- 
tion with inevitable ambiguity of the point splitting procedure. One such regularization 
is given by the lattice gauge theory, but the loss of the space symmetry makes the lattice 
loop equation both ugly and useless. Even the beautiful Eguchi-Kawai reduction could 
not save the lattice loop equation. It did not offer any alternatives to the lattice gauge 
theory. 

2.2 Momentum Loops 

Let us now forget about all the lattice artifacts and come back to continuum theory. The 
most natural thing to do with the delta function is to switch to momentum space. The 
parametric invariant functional Fourier transform of the Wilson loop is defined as follows 

W[p] = I DCW[C] exp (^-z £ dsp^{s)x'^{s)^ , (25) 

where 

DC= (f[dHs)] S'^ixiso)). (26) 

The 6 function here fixes one arbitrary point x{sq) = at the loop, the choice of Sq 
does not influence observables in virtue of translational invariance. Note, that the mass 
center fixing condition /q dsx{s) = is not parametric invariant, and would require the 
Jacobian. This is not so difficult to do, but our choice of one point fixing is simpler. 

As it was noted in[|l^] this measure factorizes for self intersecting loops on the right 
hand side of the loop equation. The argument goes as follows (assuming ti < t2, and 
choosing Sq = ti) 

(f[ dMs)] S\x{h)) (5^ (x(ti) - x{h)) = (27) 

^T<t2 \ /o-<ti+l \ 

n d'^xir) 6\x{t,)) n dM^) ^'(^(^2)) = DC,2DC2i. 

The 6 function for the loop self-intersection, 6'^ {x(ti) — ^(^2)) was combined with the 
point fixing 6 function 6'^{x{ti)) to provide the point fixing S functions for both petals. 
Note that this fits nicely large A^ factorization of the Wilson loops. For the abelian theory 
such factorization of the measure would be useless. 

The functional Laplace operator on the left hand side of the loop equation transforms 
to the momentum space as follows 

dh—-— / dh^-— => - / dhp'^ih) / dhp'^ih). (28) 

OX^{ti) Jti-O 0X^{t2J Jo ^ Jti-0 ^ 



As was discussed at length in||T2[, for continuous momentum loops this integral vanishes. 



but for the relevant loops with finite discontinuities, Apk = p{tk + 0) — p{tk — 0), the 



Laplace operator reduces to the sum of their squares J2k i^Pk) ■ These discontinuities 
are induced by emission and absorption of gluons. We shall discuss this relation in more 
detail later. 

In virtue of periodicity and symmetry this can also be written as the ordered integral 

'dhp'^ih) f"^' dt,p'^{t^) = (29) 

Jti-O 

1 rU+l 



-I dhp'^ih) J^^' dhp'^ih) 

- f^ dUp'^iU) f^^dhp'^ih) 

- I' dtip'^ih) J^' dt2p'^{h) = 
-2^ dhp'^ih) fj dUp'^iU) 



We should order the integration variables on the right as well, to avoid the trivial point 
ti = ^2 (one could trace this prescription back to the definition of the ordered exponent as 
an infinite product of the linear link factors 1 + A^(x)rfa;^, which can be varied only once 
with respect to A^{x)). After all this, we write the momentum loop equation as follows 

W[P] f dt,p'^it,) f" dhp'^iU) = f dh-^ r dh-^W [pu] W hi] . (30) 
Jo ^ Jo ^ Jo 0Pf,{t2) Jo op^[ti) 

The parts pi2,P2i of original momentum loop are parametrized by smaller proper time 
intervals T12 = ^2 — ti,T2i = 1 + ti — ^2- This can be rescaled to (0, 1) interval in virtue 
of the parametric invariance. With the manifestly invariant Anzatz (^), such rescaling is 
unnecessary. 

Note that the momentum loops, unlike the coordinate loops, are not closed. There 
are inevitable discontinuities, in particular, both parts break at the matching point. Even 
if we would have assumed the initial loop p{s) to be continuous, the W functionals on 
the right hand side would enter with the broken loops. Therefore, we have to allow such 
discontinuities on the left side as well, to get the closed equation. 

The normalization of the momentum loop functional was chosen to absorb the bare 
coupling constant. The original normalization iy[C]a;(.)=o = 1 results in the relation for 
the bare coupling constant 

Ngl = JDpW[p\ (31) 

(with the same measure Dp as for the coordinate loop C). 

We could treat this relation as a definition of the bare coupling constant. Here is where 
the ultraviolet divergences show up. We should cut these momenta, say, by the Gaussian 
cutoff exp f — e/(} (isp^(s) j. Surprisingly, such a brute force regularization preserves the 
gauge invariance, as our momenta were not shifted by a gauge field (see 1121). 
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This resembles the solution of 0{N) a model in two dimensions, where the {niUi) 
propagator was found to be the free one, with the mass given by the expectation value 
(o") of Lagrange multiplier for the constraint njn, = 1. This expectation value was then 
determined by the equation 

( UiUi ) = NgQ 



47r2 A;2 + ( a ) 

In both cases, one could regard this normalization condition as definition of unobserv- 
able bare coupling constant in terms of the ultraviolet cutoff. For the observables one 
could just leave the physical mass scale as the adjustment parameter. 

The possibility to eliminate the bare coupling constant plus the scale invariance of the 
loop equation W[p] =^ a~'^W[ap] makes this equation manifestly renormalizable. The 
physical mass scale cannot be found from the loop equation, it should rather be fitted to 
experiment. Then, just for the internal consistency check, the bare coupling constant can 
be computed as a function of the cutoff and the physical mass. 

2.3 Planar Graphs in Momentum Loop Space 

This is not to say, that the perturbation theory is unrecoverable from the momentum 
loop equation. The reader is invited to follow the direct derivation of the Faddeev- Popov 
diagrams from this equation[0, ^, starting with 



Wo[p] = Ng',l^/[p{.)-q] (32) 

where 

m-)-Q]^i[S'ipis)-q) (33) 



s=0 

is the functional 5 function. 



In higher orders there are higher functional derivatives of the functional delta functions 
s (t 'i ^Ip(-) ^ *?']■ '^^^ integration variables qi serve as the momenta inside the loops of 



planar graphs (see|]T3[ for details). 

Say, in the second order we would have two momenta gi, g2 in two windows of the one- 
gluon graph on a disk. The difference k = qi — q2 serves as the gluon momentum. There 
are two discontinuities ±k of the local momentum p{s at the edge of the disk, therefore 
the operator L{p) = 2k'^ when apphed to product VFo[pi2]M/o[p2i] on the right side. 

There are also some terms coming from the commutators of the functional derivatives 
^ with L{p). These terms drop in this order after integration over the loop. So, 



<5p(ii,2) 

when 

product VFo[pi2]M/o[P2i], it yields the gluon propagator ^. 



when inverted operator jj-^ commutes with these functional derivatives and reaches the 
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In the third order we can no longer neglect the commutator terms. These terms are 
equivalent to the [A^, A,^] terms in a field strength. They produce the triple vertex in the 
planar graphs. It is amazing, that all these planar graphs with all their ghost and gauge 
fixing terms are generated by such a simple mechanism. 

Each term is very singular, but as we argued in[0, these are spurious singularities. 



arising due to the nonperturbative origin of quark confinement. In fact, in a massive 
theory, there is no reason to expect singularities at zero momenta. The expansion coeffi- 
cients in front of the momenta grow as exp [c/qq) which is why there are singularities in 
the perturbation expansion. 

If we assume that there is an expansion in powers of momenta, the natural mass scale 
M would emerge from the zero order term W{0) = M^. The precise meaning of this 
expansion is established in the next section. We could use the scale invariance of the 
momentum loop equation to express the remaining terms in powers of W^(0). 

2.4 Concluding Remarks 

Another point to clarify. On the right hand side one should differentiate both W func- 
tional, according to the identitie^ 

^ ^ ^ ^ ^ (34) 



The resulting four terms reduce to two inequivalent ones 

2W'W' + 2W"W. (35) 

One could introduce the area derivative in momentum loop space 

^ KW^Atv, (36) 



%(t) "^"^ ' 6a^,it) 

which would make the parametric invariance explicit. 

So, finally, the momentum loop equation takes the form 

"1 rh /-i /■*2 



W[p\ 5^^ J dhp'pih) J ' dhp'^ih) = J dt2p'p{t2) J ' dtMh) (37) 



SW[pu] 5W[p2i] ^p^r . S'W[p 



'21 



Sa^aih + 0) Sa^f3{t2 + 0) ' Sa^^ih - 0)Sa^f3it2 + 0) 

We could depict it as follows 

■^Another way to arrive at the same point splitting prescription is to follow the Feynman's prescription 
of placing Afj_{x) in the middle of the interval dx in the loop product. Then dx = x{t + \dt) — x{t — ^dt) 
which leads to the half of the sum of the left and right velocities dx — ^x'{t + 0)dt + ^x'{t — 0)dt. The 
velocities x' then transform to the functional derivatives in p space. 
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+ 




It is tempting to assume that the integrands on the both sides coincide, but this is not 
necessarily true, as the perturbation theory shows. There are some terms, which vanish 
only after the closed loop integration. In the one gluon order these are the /c^/cjy terms in 
the gluon propagator. In general, these are the terms, related to the gauge fixing. 



3 Models of the Momentum Loop 

In this section we discuss various models and limits of the momentum loops to get used to 
this unusual object. In particular, we consider the QCD2, where the loops are classified 
by topology, and we also consider perturbative QCD in arbitrary dimension by WKB 
approximation. 

3.1 "Exact" Solution 



Let us first resolve the following paradox. There seems to be an "exact solution" of this 
equation of the form 

W,[p] = Zexp {^G^^jp^dp^i . (38) 
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Substituting this Anzatz into the right hand sidewe use the relations 

5W.[p\ 



G^uW,[p\ (39) 



and note that the tensor area is additive 

j Pfidpu = J Pi^dp^ + / p^rfp^, (40) 

where P12 and P21 are parts of original loop P, closed with the straight lines. These 
straight lines cancel in the tensor area integral, as they enter with opposite orientation. 

The exponentials trivially reproduce the left hand side, as for the factors in front, they 
yield 

2^2 f^ dhp'p{t2) fj dtMu)G^aG^f5. (41) 

This would agree with the left hand sideprovided 

'^ZG^a.G^P = 6a/3- (42) 

The rotations reduce the antisymmetric tensor G^j, to the canonical form 

G2k-l,2k = —G2k,2k-1 = =1= r— : ■ (43) 

with the rest of elements vanishing. 

Clearly, something went wrong. The formal Fourier transformation would bring us 
back to the constant abelian field, which cannot solve QCD. 

When we take a closer look at this Fourier transformation, we find out what went 
wrong. The Gaussian integral reduces to the solution of the classical equation 

G,Ms)=tx'^is), (44) 

which implies that p{s) could have no discontinuities, since the loop x{s) is continuous. 

In order to reproduce the gluons we must satisfy the loop equations for the discontin- 
uous momentum loops as well. In gluon perturbation theory, the Fourier integral involves 
the sum over number n of discontinuities and integrals over their locations Sj 

00 „ 
W[C] = Y, dsi... dsj (si < S2 . . . < s„) (45) 

j (/^P),,,...,,„ Ws,_sAp] exp [i jp,{t)dx,{t)) . 



These amplitudes Wn\p\ are related to the the gluon n— point functions G^^...^^ (/ci, . . . kn) 
as follows 

Ws„...,sAp\= (46) 

A A 

G^i.../.„ {ki, ...kn)S(y2 ki) « "t p^ . . . T ^^S \p{t) - J2 hO{T - s, 
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Clearly, the sum here could diverge. We expect the nonperturbative W[p] to be a smooth 
functional of momenta (nevertheless these momenta could be non-smooth functions of 
their proper time variables). 

In the case of the constant abelian field, performing the Fourier transformation from 
X to p space with n discontinuities Ap{si) = ki we would get 

/ Dx exp ( -t F^^ j> x^dx^ - i <j> Xf,dpi_, - t ^ x^(si)A;j j (47) 

oc exp (--t {F~^)ij,^ j Pf^dp^j Yl 5'^{ki). 

The product of the delta functions comes from the integrations over the positions of Xi. 
This singularity arises because we have no iP' term in the exponential to damp spikes in 
x{s). 

So, the simple exponential in momentum space does not correspond to the constant 
abelian field. It would rather correspond to violation of the gauge invariance, as we could 
only obtain it by integrating over the discontinuous loops in coordinate space. Therefore 
we have to reject this formal solution. 

3.2 Perturbative QCD 

Let us now discuss the perturbative solution. The formal perturbation expansion for 
the momentum Wilson loop does not reflect its asymptotic behavior. To get a better 
description of the asymptotic freedom in momentum loop space, we should rather expand 
in the exponential of the path integral 

W\p\ oc / Dxexp {\og{yV{x)) - I i p'^{s)x^{s)\ , (48) 

near some (imaginary) saddle point x{s) = C{s). It should satisfy the equation 

^Pli^) = CUs)'-^f^. (49) 

Ocr^^[s) 
In the first nontrivial order we get an equation of the following form 

ip'^.{s) = Ng\C) Cl{s) j {dC,{t)d, - dC,{t)d^) D (C(s) - C{t)) . (50) 

Here the running coupling constant Ng'^{C) corresponds to the scale of C and 

is the gluon propagator. We could take this classical equation as a parametrization of 
Pc{s)- We observe that pc ~ Ng'^{C)\C\^^'^. In four dimensions, for the small loop C, 
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the corresponding momenta are large. In two and three dimensions, large coordinate loops 
correspond to large momentum loops. 

As for the powerlike divergency of the contour integral at t ^ s in three and four 
dimensions, it could be regularized by the principle value prescription x^ ^> x^ + e. For 
the smooth nonintersecting loop C such integrals are finite. 

The t § C^dpQ term reduces to —2 times the same one gluon exchange graph in virtue 
of the classical equation. So, we get the following large p asymptotics for W[p] in the 
C-parametrization 

\ogW[pc] = +\Ng\C)jdC,{s)jdC,{t)D{C{s)-C{t)) ~ Ng\C)\C\'~'. (52) 

Unlike the formal perturbation theory, this expression is a well defined parametric invari- 
ant functional of C{s). One could go further and compute the quadratic form for the 
second variations 5C{s) near the saddle point. This would give the pre-exponential factor 

inW[Pc]. 

Note that in two dimensions there is no pre-exponential part for the non-intersecting 
loop. The first order of perturbation theory is exact for logW^(C) in this case. The 
integrals can be done, and they reduce to the (absolute value of) the tensor area. Thus, 

p''c{s)=iNg\,,C''{s), (53) 

and 



logW^[P] 



j Pidpi 



2Ng^ 
up to the nonlocal corrections due to self-intersecting C-loops in the path integral. 



(54) 



At present we do not know how to compute these corrections, but at large pc the 
fluctuations 5C ~ } are small compared to the classical loop |C| ~ -^^. Therefore, 

for large p we have the area law for W^[p]. This is the general phenomenon, as one can 
see from the momentum loop equation. In fact, the above "derivation" of the tensor area 
law is valid for the large loops, where discontinuities and self-intersections are irrelevant. 



3.3 Beyond Perturbation Theory: Random Vacuum Field 

To get a rough idea what could be the momentum Wilson loop beyond perturbation 
theory, let us consider a simple "model" of QCD vacuum as random constant abelian 
field F^i, with some probability distribution P{F). As we discussed above, this is not a 
solution of the loop equations. Such model misses the hard gluon phenomena, responsible 
for asymptotic freedom, but it is the simplest model with confinement. 

The Wilson loop would be 

W{C) ~ J P{F)dF exp i^i F^, j x^ dx}j , (55) 

16 



and the corresponding momentum loop (for continuous p{s)) 

W{p) ~ J P{F)dFVdetF exp (^-^i F^J j> p^ dpj^ . (56) 

The pfaffian -y/det F comes as a result of computation of the functional determinant by 
the ( regularization in Appendix A. 

Both loop functionals could be expanded in Taylor series in this simple model. The 
Taylor expansion at small X exists, because there are no gluons, just the constant field. 
This is a wrong feature, to be modified. The Taylor expansion at small p exists, because 
the Wilson loop decreases fast enough in this model. This feature we expect to be present 
in the full nonperturbative QCD. 

Note, that the Fourier integrals here are not the usual integrals with the Wiener 
measure, corresponding to Brownian motion. These are the parametric invariant integrals, 
without any cutoff for higher harmonics of the field x{s) or p{s). Such integrals can be 
defined as limits of the usual integrals, with discrete Fourier expansion at the parametric 
circle. 

oo 

x(t) = ^a„e2-"*. (57) 

oo 

For the readers convenience, we consider these integrals in Appendix. We also compute 
the divergent pre-exponential factors with the ( regularization. 

The important thing to know about these integrals is that they exist. One does 
not have to break the parametric invariance by first adding the gauge fixing terms and 
then compensating that by some Jacobians. One could simply compute them with some 
invariant regularization, such as (, or Pauli-Villars regularization. 

This also regularizes the wiggles, or backtracking. The Wilson loop does not change 
when a little spike with zero area (backtracking) is added to the loop. Such backtracking 
are hidden in the high harmonics a„ of the position field. The regularization of these 
wiggles leads to finite result in the constant field model. We expect the same in full 
nonperturbative QCD. 

In this case, the position field could exist, unlike the usual master field. Our Anzatz 
(I) yields for the Taylor expansion terms 

Wn{p) « ;^trt (^f^dp,is)X,y = (58) 

j^ j ■■■ j d{si<S2...<Sn) dp^^{si) . . . dp^„{Sn) tr (X^, . . . X^J . 

These terms, in fact, reduce to traces of products of commutators, as there is translational 
symmetry X^ =^ X^ + 1 * a^. 

Comparing the above path integral for these terms with the second quantized version 
we conclude that the string correlators 

^J Dx6\x{s^)) W{C) x^,{si) . . .x^„(s„) (59) 
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do not depend of proper times Sj, as long as these times are ordered. In general case 
these are combinations of the step functions, as they should be in parametric invariant 
(i.e. ID topological) theory. Integrating these terms by parts with dp^j^{si) . . . dp^^{sn) 
we would reproduce the commutator terms coming from derivatives of the step function 
6 {si < S2 ■ ■ ■ < Sn) in the ordered operator product representation. 

These string correlations correspond to the Neumann boundary conditions. Here we 
do not fix the edge C of the " random surface" , but rather let it fluctuate freely. This fixes 
the corresponding local momentum: p^ = at the edge of the surface. 

Adding the ^ax'^ term at the edge would correspond to the term |^p^ in the Hamilto- 
nian. We do have such terms in our physical Hamiltonian, of course, but for the purposes 
of computation of the simplest correlators we could set a — » 0, which is equivalent to the 
boundary condition p = 0. 

Our approach to the "QCD string" is similar to the matrix model approach to the 
2D quantum gravity. We shall work directly in parametric invariant "sector" rather than 
fixing the gauge and then imposing constraints to come back to the physical sector. 



4 Equations for Dual Amplitudes 

The momentum loop equation is still ill defined. The meaning of the functional derivatives 
j-Tj-T at the breaking points ti,t2 needs to be clarified. Also, the necessity to include 
momentum discontinuities makes the discrete Fourier expansion useless. Ideal thing would 
be some nonsingular equation with functions of finite number of variables instead of 
functionals. 



4.1 Kinematics 

Let us try to obtain the closed equations for the dual amplitudes 

A{h,...K) = j^tTV{h)...V{K), (60) 

where 

V{k)=exp{tk^X^) (61) 

is the vertex operator, and momenta ki add up to zero. 

These traces are cyclic symmetric, of course. In addition there are relations between 
amplitudes with various number of points 

Ai...,^k,vk,...)=Ai...,{^ + r])k,...) (62) 

which do not have analogies in the particle theory. In particular, the two point function 
is trivial 

A{k, -k) = A() = 1. (63) 
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The first nontrivial one is the three point function A(k, q — k, —q). Given the generic n 
point function, we automatically get all the previous functions by setting some momenta 
to zero. 

The momentum loop represents the polygon with the sides ki and vertices 

i 

Pi = Yl ^v ki=Pi- Pi-i,Po = 0. (64) 

In parameter space the sides are the discontinuities, they take infinitesimal time steps. 
Si < t < Si + e. This does not matter, in virtue of parametric invariance. The only thing 
which matters is the geometric shape of this polygon. 

These amplitudes corresponds to W[p] in the limit 

n 
pis)^Y.k^Ois-S^), (65) 

when there is no smooth part of momentum. Such limit is quite natural from the point 
of view of the string theory, but in perturbative QCD it would be singular. 

Existence of such limit implies convergence of the path integrals with W{C), which 
is violated in formal perturbation theory. In above constant field model such limit exist. 
The symplectic form becomes the quadratic form in momenta 

/n i—1 

P,dp.^Y.F;uKllk^ (66) 

Given the probability distribution function for the tensor F one could compute the cor- 
responding amplitudes in this model. 



4.2 Loop Equations 

Let us turn to the loop equations in real QCD. We already know the limit of the left hand 
sideof the momentum loop equation. It is just 

L{W){p) -^ (jZki^ A{k„...k^). (67) 

This involves the sum of squares of lengths of the polygon. 

The problem it to find the limit of the right hand sideof the loop equation. We have 
to compute it for the general Anzatz(^), and then tend the local momentum to the sum 
of the theta functions. 

The double integral around our polygon breaks two of its sides 
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The extra discontinuity of momentum 

A=p(ti)-p(t2) (68) 

arising with opposite signs in the two momentum loops leads to extra vertex operators 

W\vr2\ = ^tr (f exp {% £' dv,X^ V(A)) , (69) 

W\v2x\ = ^tr (f exp [i £' rfp^X^) \/(-A)) . 

These extra vertex operators correspond to internal sides, closing these two loops. They 
are not parts of external loop p, so they are not involved in the functional derivatives. 

Let us now apply the functional derivative ^j-jv to our Anzatz. Start with the identity 
5 



— I ■ 



Texp U / dp^xA = Texp U / dp^xA X^^Texp U / dp^xA , (70) 



6p'^{u 

and differentiate it in u. We find 
6 



Spuiu) 



Texp U / dpaX 



t-:^ ^ . , . T exp ( 2 / dp^X^ 



du 5p'j^{u) 



(71) 



p'^ (u) fexpit dp^xA [Xy, Xa\ f exp 



dp^X^ 



These are momentum space analogs of the Mandelstam relation for the area derivative, 
with X^ playing the role of the gauge field. These identities play the crucial role in our 
construction. Without them we would have only functional equations, for the general 
function p{s). 

When the momentum tends to the sum of the theta functions, the integrals over ti, ^2 
are concentrated inside the infinitesimal intervals (sj, Sj + e), (sj, Sj + e) , where derivative 
p' is present. Due to our principle value prescription they can never fall in the same 
interval. Inside such intervals the functional derivatives split the vertex operator 

5V{h) oc r^^dp^it) (\/(A;0[X„Xj ® V{k':) + V{k'^ ® [X„Xj\/(A;f)) . (72) 

J Si 
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where 

K = Pit) - Pi-s^), k': = Pis, + e) - Pit), k[ + k'l = h. (73) 

The first factor in the tensor product goes in one trace, the second one goes in another 
trace, since the loops reconnect precisely at the variation point t = ti or t = t2- The sum 
of two terms with commutators switching from one trace to another arises because our 
definition of velocity x'iu) involves the sum of left and right derivatives. 

Inside each infinitesimal interval dpit) = kd^, k' = C,k, k" = (1 — ^)k, where ^ = i^^^. 
Therefore we get 

5Vih) oc k^£ di {yiih)[X,,X^] ® \/((l - i)ki) + Viih) ® [X„XJ\/((1 - i)h)) . 

(74) 
The new objects here are the commutators. However, they can be obtained as a limit of 
the vertex operators 

[X^, XJ k^ = 2 (^^Vipk)Viq - pk)Vi-q)\ (75) 

Geometrically this corresponds to adding little triangle in the /x direction as shown at the 
above figure. 

Now we have all the ingredients in place. We obtain the following equation 



A ir]kj, kj+i, . . . , (1 - ^ + pi)ki, gi - pi/cj, -gi. A) 

A i^ki, ki+i, ...,il-r] + p2)kj, q2 - P2kj, -q2, -A) + 

A (-gi, gi - pih, (^ + pi)ki, ki+i, . . . , (1 - r] + p2)kj, g2 - p2kj, -ga, -A) 

Air]kj,kj+i,...,il -^)ki,A) 

= 0. 



(76) 



where 



A = efc^ + . . . + (1 - v)k, = - ivk, + . . . + (1 - Oh) (77) 

is the transfer momentum. 

These are exact equations for functions of finite number of variables rather than func- 
tionals. The main advantage of this equation over the old loop equation is absence of 
singularities. There are no delta functions here, the only /c-dependence is polynomial. 

5 Taylor Expansion at Small Momenta 

One may expand these amplitudes in momenta (which should be possible in massive 
theory), and recurrently find the expansion coefficients. An important property of these 
equations is the fact that the same coefficients serve all the amplitudes, for arbitrary n. 
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5.1 Equations for Commutators 

These expansion coefficients are related to the traces of products of the X operators. 
Such traces are universal tensors, composed from 5^^, independently of the number of ki 
variables multiplying these tensors. 

We shall denote by Al the O (k^j terms in the amplitude, for generic momenta 
ki, . . .kn- In fact, this Al would involve only combinations of L or less different momenta, 
so that remaining momenta are redundant. This allows us to find relations for these 
coefficients, taking only finite number of momenta. 

The leading terms are quadratic in momenta. Matching these two quadratic forms 
gives the ffist nontrivial relation. The A term yields J2 kf- In the remaining A^ A terms 
only the product of the zeroth and fourth order terms A4 contribute. The fourth order 
terms reduce to the trace of two commutators 

— tr [X^, Xa] [X^, Xp] = a (S^J^p - S^pS^a) ■ (78) 

We find the following sum ruleQ 

2{d-l)a = l (79) 

In higher orders we get more relations, each time involving the higher derivative terms, 
i.e. higher commutators of position operator. In order to study these relations systemat- 
ically, let us parametrize the conserved momenta as ki = pi — pj_i and rewrite the dual 
amplitude as follows 

^(Pl -Pn,---,Pn-Pn-l) = ^tr {S{pu Pn) ■ ■ ■ S{pn, Pn-l)) (80) 



with 



S{k,q) = V{q)V{k-q)V{-k)= f exp (t J^X^dp^^ (81) 



where the ordered integral goes around the triangle made from three vectors q,k ~ q, —k. 
These triangles ffil the polygon with vertices Pi and sides k^. 

At small momenta this formula reduces to 

Sik, q)^l + h^q, [X^, X,] + '- (k^k.qx - q^qukx) [X^, [X„ Xx]] + ■ ■ ■ ■ (82) 

This is the nonabelian version of the Stokes theorem for the triangle. One can find higher 
terms here using the ordered integrals. From this representation it is obvious that the 
Taylor expansion starts from the quartic term, as the trace of commutator vanishes. 

In higher orders we get the traces of polynomials to the set of traces of all possible 
commutators 

Xu^u2...vn = [Xui, [Xu2, ■ ■ ■ X^^] . . .]. (83) 

"^It was first obtained by different method in collaboration with Ivan Kostov(unpublished). 

22 



The general translation invariant trace would involve products 

These T-traces are not independent tensors. Those would be given by the usual traces 

G(/ii.../i„) = — trX^, ...X^„, (85) 

except these traces do not exist because they are not invariant under translations 6X = 
a*l. 



5.2 Planar Connected Moments 

The following parametrization proved to be convenient. Let us introduce the planar 
connected moments VF^i...^„ defined recurrently as follows 

W, = Gifiy, (86) 

W^Mi-Mn =G{fii...fin)- E ^m(i) • • • ^/^(-f • 

planar 

Here the sum goes over all planar decompositions of the set (/xi . . . /i„) into subsets 
f/i'^-'^M , . . . (/i*^'"M. For example, 

W^, = G{fiu)-G{fi)G{u), (87) 

W^^x = GifiuX) - G{fi)G{u\) - G{u)G{\fi) - G{\)G{fiu) + 2G{fi)G{u)G{\). (88) 

In the next section we shall use this representation to build the Fock space representation 
of the X operator. 

The advantage of these C— moments is their simple behaviour under translations. 
Only the first one gets translated 

6W^ = a^ (89) 

while all the higher ones stay unchanged. This can be proven by induction. The variation 
of the G— moments under infinitesimal translation reads 

5 G {ni . . . fin) = Y^ a^, G (fii+i . . . in^i) , (90) 

I 

with the cyclic ordering implied. Straightforward planar graph counting gives the same 
for the variation of the planar representation 

G (/ii . . . /i„) = W,,,„,„ + Y. W^M(i) • • • ^m(™) ' (91) 

planar 

5G{fli...fln)= J2 ^I^W ■ ■ ■ ^W^W^/.(™) = 
planar 

I planar 



/ , a,^^G {fii+i ■ ■ ■ fii-i) 
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In the next section we present a simpler proof based on operator representation. 
So, we could fix the gauge by imposing 

W^ = (92) 

and parametrize the traces by the invariant coefficients in front of various terms in the 
planar connected moments. Those are constant tensors, composed from 5^-^^ . To simplify 
the formulas we shall use the notation 

8i,j = 5^^^^^. (93) 

This shall never lead to confusion, as all our indexes would be summed over at the end 
anyway. 

Only cyclic symmetric tensors may appear. Up to the eighth order 

W^i,2 = h,2 A, (94) 

W^l,2,3,4 = <^1,3 ^2,4 B + ((5i,4 ^2,3 + ^1,2 <^3,4) C, 

W^l,2,3,4,5,6 = (^1,4 ^2,5 hfi D + 

((^1,4 (^2,6 (^3,5 + (^1,5 ^2,4 (^3,6 + (^1,3 (^2,5 ^Afi) E + 

((^1,6 (^2,5 ^3,4 + h,2 (^3,6 (^4,5 + (^1,4 h;i ^5,6) F + 

(5l,5'^2,6'^3,4 + '5l,6'^2,4'^3,5 + <5l,3'^2,6'^4,5 + '^1,5'^2,3^4,6 + '^1,2^3,5'^4,6 + '^1,3^2,4'^5,6) G + 

(^1,6 ^2,3 (^4,5 + '^1,2 (^3,4 ^5,6) H. 

From the above solution for the fourth moment we find 

^ = ^ + '^'+2(d^- (^^) 

The coefficients A, C remain arbitrary. 



5.3 Higher Orders 

In higher orders, let us pick out the highest O [k^j order term in (|76|), moving the lower 
ones to the right 

Y, kt k^ f di f d^ (96) 

^^. Jo Jo 

ti{[X,,X^]V{A)[X,,Xp]Vivkj)Vik,+,)...Viil~Oki))^,_, 
= F{Ai<L_2,kj). 

The left side involves (linearly!) the highest C coefficients, with L + 2 indexes. The order 
of the coefficients on the right side is L and lower. 
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So, this chain of equations have triangular form, which is a great simphfication. The 
number of equations we get in the O (k^) order is equal to the number of the L-th order 
cyclic symmetric invariants made from n momenta. The number of unknowns coincides 
with the number of cyclic symmetric tensors with L + 2indexes, made from 6^^. 

There are, in general, much more equations than unknowns, but some hidden symme- 
try which reduces this system to just one equation, at least up to the orders we computed. 

Thus, all the O {k"^) equations for any number n of legs reduce to (PBI). The O (fc^) 



equations we studied, with up to 6 legs, all reduce to the following equation 

3D = 2A^ + A l—^^—-4:C]+E-4:F + 7G-H. (97) 

\2{d-l) ) 

These equations were obtained by means of Mathematica . We generated all cyclic sym- 
metric tensor structures recurrently, then introduced coefficient 

in front and substituted into the operator expansion. The Cuntz algebra was imple- 
mented to compute the expectation values of the products of the X operators. 

We had to use some Mathematica tricks to increase efficiency in computations of the 
loop equations. In its present form the program could compute a couple more orders with- 
out much effort at the personal workstation. The program is included in the uuencoded 
hep-th files and could also be provided by e-mail upon request. 

In the next order we found the equation, but it is too lengthy to write it down here. 
It is included in the Mathematica file. 

We could not find the general formula for the number of parameters in d dimensions. 
The group theory must have these answers. Clearly, one cannot alternate more than d 
indexes, which restricts invariants made from d + 1 or more 5^i, tensors. 

In four dimensions this would start from ten indexes which we do not consider here. 
However, in two dimensions these restrictions are quite important. 

5.4 Symplectic Symmetry in QCD2 

There is, in fact more symmetry in two dimensions. The symplectic invariance of 2D 
QCD corresponds to the following transformations 

Sx^{t) = ae^^Xy{t), (98) 

^P^l{^) = -ae^^p^{t) 

One can readily check that these transformations preserve the symplectic form 



j> Pfidx^ 



as well as the tensor areas 



f e^:,Xf,dx^, j e^„Pf,dp„. 
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We already saw this symmetry in the WKB approximation at large p when the momentum 
loop was a function of the tensor area in momentum space. 

These transformations lead to the symmetry of our equations in 2D 

6X^ = ae^,X,. (99) 

As a consequence of this symmetry, one should use only the e^j^ symbols in constructing 
the expectation values of products of commutators of X. The parity implies that the e 
symbols come in pairs, therefore only the multiples of four X operators could have nonzero 
traces. The lowest commutator X^^ = [X^,X;^] is symplectic invariant, but the higher 
ones, in general, are not. This leads to extra restrictions on the invariant structures, to 
be used in the solution. 

Up to the eighth order the symplectic invariant solution in two dimensions reads 

A = 0, B = l, C = -]:... (100) 

6 

It corresponds to 

5.5 Higher Dimensions 

In higher dimensions it is tempting to try something similar. Due to the Z4 symmetry 

p^{s) ^Ljpf,{s), u* = 1, (102) 

of the momentum loop equation there is always a solution with Vri,...^„ = unless n = 
OmodA, so that the Taylor expansion goes in p^. 

This is not sufficient to eliminate the higher invariants in C— tensors. Say, in the 
fourth order the second invariant C would still be present 

Similar terms would exist in higher orders. We do not have symplectic symmetry here to 
reduce invariants, so what do we do? 

In principle, we should fix the unknown coefficients by matching with perturbative 
QCD at large momenta. This is easier said than done. As a model, we could try to build 
the general terms from commutators X^^^ in arbitrary dimension. Say, in the fourth order 
this would correspond to 

C = -\b, (104) 

so that the C— tensor would be equivalent to the double commutator 

W^l,2,3,4 = ^ -B (2 5i,3 ^2,4 - ^1,4 ^2,3 - h,2 ^a) ■ (105) 

26 



Solving the loop equation we find 



^ = WTTTy '^°^' 



In four dimensions we could get some number of Taylor terms analytically. Then, 
the Pade approximation could be used to estimate the spectrum. One could also try to 
compare some string models with these relations. 

In Appendix B we discuss the possible numerical solution of these equations. We post- 
pone the numerics till the better understanding of analytic properties of the momentum 
loop equation. 

6 Operator Representation 

The equations of the previous section can be used to recurrently determine the planar 
connected moments PV^^...^^. 

Suppose we solved this algebraic problem. The natural next question is how to recon- 
struct the matrix X from its moments. In the large N limit this problem can be explicitly 
solved by means of the noncommutative probability theory. 

6.1 Noncommutative Probability Theory 

The noncommutative probability theory in its initial formj^ covered only so called free 
noncommutative variables, corresponding to the factorized measure. The theory of inter- 
acting noncommutative variables which we use here, was developed in the recent work|TT|, 



where it was applied to the matrix models and QCD2. Some elements of this theory were 
discovered back in the eighties in the context of the planar graph Schwinger- Dyson 
equations. 

Let us introduce auxiliary operators a^ and aj^ satisfying the Cuntz algebra 

Qf^al = 6^u (107) 

which is a limit of the q— deformed oscillator algebra 

afMal-qala^ = 5^^ (108) 

when g -^ 0. These operators correspond to the source j = a^ and noncommutative 
derivative j^ = a introduced by Cvitanovic et. al.p). 

We could explicitly build these operators as the limit of the large N matrix aj^ and its 
matrix derivative a^ = j^jtt- The trick is that we shall keep all the traces of products of 

a^ finite in the large A^ limit, so that the moments vanish except the zeroth 

l^tTal...al^=6no. (109) 
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The matrix derivative is the N = oo counterpart of the matrix of the derivatives with 
respect to the matrix elements of a^ 

The transposition of indexes ij is needed to obtain the simple matrix multiplication rules 

(a.al ...al) = {aS' (al . . . aH = (111) 

^uui ( c^L • • • ctt ) + TT y^ ^uu, tr ( a|„ . . . a|, ) ( a|, . . . a|, ) -^ 5uu, ( aL • • • a|, ) 

In the first term the factor of A^ in definition of the matrix derivative cancelled with the 
trace ba = N. The remaining terms all contain the vanishing moments. 

6.2 Operator Expansion for the Position Operator 

Let us define the following operator expansion 

X^ = W,*l + a, + Y: W,,„„,^al ... 4^ (112) 

where the sum goes over all planar connected moments C with all tensor indexes. The 
first term is a c-number. It takes care of the possible vacuum average of the position 
operator. It drops in the dual amplitudes, as it commutes with the rest of terms and 
corresponding exponent can be factored out of the vertex operator 

V{k)=exp{tk^X^) -^exp{tk^W^)V{k). (113) 

These phase factors all cancel due to momentum conservation, which allows us to drop 
the W^— term from the position operator. 

With this representation it is obvious that the planar connected moments W,,, are 
all translational invariant, except for the first one, Wfj,. The translation of the c-number 
term provides the correct shift of the position operator, and there is no way the variation 
of higher terms could have compensated each other, as they enter in front of different 
products of the aj^ operators. 

Consider the expectation values of the products of such operators with respect to the 
"ground state" |0) annihilated by a^ 

ajO) = 0, (0|at =0. (114) 

With our matrix derivative representation of these operators 

|0) = 1, (0| = (5(at); (115) 
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so that 

(0| a,,... a,„F{a^) |0) = 1-^ . . . ^.^F{a^ ^ 0) (116) 

Clearly, only the terms with all a and a^ mutually cancelled will remain. This gives 
the usual decomposition of the moments into planar connected ones 

(0| X,,... X,, |0) = W,,...,^ + Y: W.W... + ... (117) 

The details of this beautiful construction for matrix models and 2d QCD can be found in 



m 



In our case, the coefficients of the operator expansion are governed by the loop equa- 
tions. Up to the sixth order 

X^ = a^ + Aal + Baialal + C (at aj^a^ + ^ aj^aj) + . . . (118) 

with the coefficients A,B,C,... satisfying above equations (p5|),(p7D. 
The simplest solution, discussed in the previous section, reads, 

1 



'^'^-"^+6(rf 



a], 



a^ a^ 



+ 0{{a^y) (119) 



What could be the use of this expansion? Perhaps, some truncation procedure using 
some trial functions can be found. Or else, one can always use the analytic methods, by 
extrapolation of expansion for the physical observables by the Pade technique. We do not 
know all the expansion parameters, they are to be fixed by matching with perturbative 
QCD at large momenta. 

6.3 Perturbative vs Nonperturbative Fock Space 

It is instructive to compare this construction with that of [O] . We have here (at arbitrary 



dimension of space) much smaller Fock space than the one found by Gopakumar and 
Gross for QCD2. There, the operators a, a^ had the momenta as their arguments, which 
made the Fock space similar to the space of the free gluons (except for the Boltzmann 
statistics instead of the Bose statistics of the usual free gluon theory). 

The present construction has much less degrees of freedom, in fact, no continuous 
degrees of freedom at all. We simply have d pairs a^, a^ which is like the Fock space of a 
single Boltzmann oscillator. The states are spanned by all d— letter words 

|/ii/i2 . . . /^n) = al^al^ ...al^ |0) . (120) 

This tremendous reduction of degrees of freedom (very much like the one in Eguchi-Kawai 
model, but in continuum theory) takes place due to implied confinement. 
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It is due to confinement that we are able to expand amplitudes in local momenta, so 
that only the Taylor coefficients need to be described. The Taylor expansion is much sim- 
pler that the Fourier expansion, where the basic set of functions are labelled by continuous 
momenta. 

Another simplification came from unbroken parametric invariance. The general Taylor 
expansion of the loop functional would involve unknown coefficient functions, which would 
result in the loop operators a^{s), like those of the conventional string theory. Expanding 
those in discrete Fourier expansion on a unit circle we would get the usual infinite set of 
the string oscillators a^,„. 

So, why is our Fock space smaller that that of the string theory? This is because we 
did not fix the parametric gauge either. Should we break the parametric invariance by 
adding the conventional x't term in the exponential in the definition of the path integral, 
(i.e. use the Wiener measure instead of the fiat one we use here), we would have all these 
problems back again. 

However, as we learned from the recent success of the matrix models, the quantum 
geometry problems can be solved by purely algebraic and combinatorial methods, without 
any breaking of parametric invariance. The (perturbative) Fock space of the Liouville 
theory in conformal gauge involves continuum fields, but we know that there is a different 
picture, with only discrete degrees of freedom. 

This picture is the matrix model quantized by means of the noncommutative prob- 
ability theory. The Fock space in this picture involves only one Cuntz pair a, aK The 
perturbative Fock space was too big as it turned out after exact solution. After that, the 
Liouville theory Fock space was also reduced, by methods of the topological field theory. 

So, we expect the same thing in the large A^ QCD. The nonperturbative Fock space 
must be much smaller than the perturbative one due to the color confinement. The 
operators we build here are discrete because we are working in purported confining phase 
of QCD, where the motion is finite, so that most of the perturbative gluon degrees of 
freedom are eliminated. 

By the same reason we switched from the original gauge field (the Witten's master 
field) to the conjugate coordinate operator. The master field was unbounded even in 
confining phase due to the hard gluons. The coordinate field is bounded in confining 
phase, which is infinitely simpler to describe in our picture than the free motion. 

In fact, without confinement, our methods would be completely pointless. The de- 
scription of the plane wave of unconfined gluon in terms of Taylor coefficients would be 
very inefficient. 

The proof of quark confinement within such approach cannot be achieved, all we can 
say is that exact equations for the dual amplitudes allow for a systematic expansion in 
local momenta, which can be interpreted in operator terms. 

Note however, that these same equations, being iterated from perturbative vacuum 
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instead, generate all the planar graphs of perturbative QCD. The equations being non- 
linear, we cannot say whether our Taylor expansion describes the same solution, which 
matches perturbative QCD. 

7 Hadron Spectrum 

In this section we derive the wave equations for the mesons and glueballs. These equations 
correspond to a different physics: stringy states described by a long timelike Wilson loop 
in meson case, and small fluctuations of the Wilson loop in the glueball case. Still, in 
both cases, the wave equations can be explicitly written down in terms of the position 
operator. 

7.1 Mesons 

The equations of the previous section are supposed to fix X up to the overall gauge 
transformation X^ =^ S^^X^S. Let us assume we found the solution for X, numerically 
or analytically. What can be said about the particle spectrum? Let us study the meson 
sector first. 

Consider the limit, when the loop becomes an infinite strip. This corresponds to 
the propagation of the meson made from quark and antiquaqrk. There are two position 
operators X^, X^, and the corresponding loop function can be regarded as tensor product 
of the two ordered exponentials 

* (r. X, i) o. / Dp(.) JDpOt exp (£ M (, X„^ - i/) ® (> i, * - h) ) 

(121) 
with the above Dirac Hamiltonians ([T0|). 

Let us shift the time T ^ T+dT and try to get the Schrodinger equation (in imaginary 
time). Introducing the vectors g^ = p^(T + dT) — Pfj,(T) we get the vertex operators 

Texp U j^ X^dp^{t)\ -^ exp [iq^X^) = V{q). (122) 

Now, expanding the rest of the factors in dT 



exp {-dT {m + i j^.Pf.iT + dT))) = exp i^-dT i^H + 1 7^'?^)) -^ I - dT i^H + 1 7^g^^ 

(123) 
and comparing the linear terms, we get the following 

- ^ (t, X, 1) = ((h + j,A,) ®{h + %A,)) ^ (t, X, 1) , (124) 

where 

d'^q 

("2^ 



4 = ^ I 7^,Vi^)l^^ = 9A^)- (125) 
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We also get some terms with the quark mass corrections, which we absorb into the bare 
mass m. 

The A^ term is the net contribution of the "QCD string". This is the natural general- 
ization of the derivative operator in the Dirac equation. Incidentally, this is what we get 
for the Witten's master field. 

For a finite matrix X this term would be a sum of the (gradients of the) delta func- 
tions, whereas in the large N limit it is related to the density of eigenvalues at small 
X. The singularities at small loops would correspond to the divergent integrals at small 
eigenvalues. 

There are d different sets of eigenvalues /t^ and the same number of different eigenvec- 
tors flfi, as the matrices X^ = fi^/t^fij^ do not commute. Therefore, the angular variables, 
corresponding to the eigenvectors of each X^ would not decouple here. 

With the operator language described in the previous section, the gradient of the delta 
function can be defined by its matrix elements with some trial functions in the Fock space 
of the operators aj^. 

Once we know the master field, the above Schrodinger equation is quite explicit. The 
states with fixed total momentum 

plot ^ pstring ^ p ^ p (126) 

are selected by condition 

tP';''''''^{T,X,]() = L (^) +tr (^]U(t,X,^), (127) 

corresponding to total translation. The mass spectrum is given by stationary solutions, 
i.e. 

((iJ + 7^4)©(i? + 7^4)) ^(X,1)=0. (128) 

The appropriate trial states are superpositions of the vertex operators 

^fc„....„,S„...fc^ {X, ^)=V ih) ...V (kn) ® V {h) ...V {k^) , (129) 

The above wave operator acts linearly on the space of such states. Note, that we are not 
considering the dual amplitudes now. The quark degrees of freedom are present in this 
hamiltonian. Still, the vertex operators can be used as the basis for the wave equation. 
The computational aspects are beyond our present scope. 

7.2 Glueballs 

In the glueball sector the situation is different. The glueball wave function G{C) can 
be obtained of as variation of the Wilson loop G{C) = 6W{C), in the same way, as the 
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variations of the Higgs condensate describe the physical excitations in the phase with 
"broken" gauge symmetry. 

This relation between G{C) and W{C) can be formally derived from the loop equation 
for the generalized loop average 

Wq{C) = (^QtT (^f exp [^j^A^dx^^ ^ (130) 

where Q is some invariant operator with the fixed total momentum, 

Q = y"d'^l/exp(-z%)tr (P(F^,,V,F^^,...)) (131) 

One could think of this Q as an infinitesimal loop, moved around with the plane wave 
weight exp {—iky). 

The functional Wq{C) is not translation invariant, it is rather the momentum eigen- 
state 

jdsj^WQ{C)=tk,WQ{C). (132) 

as it follows from the translation y ^ y + a in above integral. 

The Levy operator still yields the Yang- Mills equation inside the trace. In virtue of the 
Schwinger-Dyson equations this leads to the usual loop splitting terms plus the contact 
terms, which we ignore, as we are interested in the residue in the pole at the glueball 
spectrum fc^ = —w?. The loop equation for this residue Gq{C) reads 

L{Gq){C) = R{Gq,W){C) + R{W,Gq){C) (133) 

We used here the factorization theorem, keeping in mind that for finite momentum k 
the expectation value of Q vanishes. Hence, the leading term in the average product of 
the three invariant operators reads 



Qtr I Texp U^ A^dxA 1 tr I Texp U^ A^dxA 1 ) -^ (134) 

Qtrftexpf£ A^dxA\\Uiifey.^U A^dxA\\ + 
Q tr T exp i A^dXf, ) ( tr T exp f A^dx^ ) 



Now, the loop equation ( |133| ) would be satisfied for the linear variation 

Gq{C) = 6W{C, X) = w[G,X + 6X) - W{G, X), (135) 

along the nonlinear loop equation. In other words, if the nonlinear loop equation is 
satisfied, when the position field X is shifted by 6X, then the variation of W{G,X) 
satisfies the linearized loop equation (|133|) . 
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Such perturbation corresponds to the variation of the vertex operator 

6V{k)= I d\V{\k)ik^5X^V{{l-\)k). (136) 

with 5X^ represented by the operator expansion 

5X, = 5W, + Y. ^W,,,...,^al^ ...al^, (137) 

The set of equations for the perturbed amphtudes 

6A{k,, ...k^) = j^ik^ f d\ tr 6X^ V{{1- \)ki) V{ki^r) . . . V{h_^) V{\ki) (138) 

would now provide the linear relations for the coefficients 5W,,,. This variation should be 
the momentum eigenstate, i.e. 

5A{ki,...kn)(x5'^{k-Y,h) (139) 

after which we get linear homogeneous system of equations, involving A;^ as a parameter. 
The solvability conditions should provide the spectral equation for Q. In virtue of the 
space symmetry this would quantize the square fc^ = — M^. 

8 Conclusion 

This is as close as we could get to the QCD string. We built the model of d large A^ 
matrices X^ representing position operators of the open string endpoint. The infinite 
set of algebraic relations for the matrix elements of of X was obtained. The Fock space 
representation of these infinite matrices in terms of the operator expansion was found. 

The expansion coefficients are recurrently related to each other, which allows them 
to be explicitly computed in QCD2, and strongly reduced in general case. The problem 
of elimination of unknown parameters in the Taylor expansion in general case remains 
opened. 

There is nothing like the functional integral over random surfaces, which is not so bad 
after all. Such integrals have problems even in critical dimension of space, where there 
are tachyons. As for the noncritical strings, they do not yet exist. Hopefully we would 
see such string theories in future ITBI, then they would have a chance to be equivalent to 



QCD. 

The advantage of the present "string" theory is that it is manifestly parametric in- 
variant. We do not reduce the larger state to the parametric invariant sector by Virasoro 
constraints, we directly work in physical sector. In this sense we are building our "string 
model" as one dimensional topological field theory. 
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We bypassed the hard unsolved problem of computation of the Wilson loops in coor- 
dinate space. Not only these loops are hard to compute, they are not defined as numbers 
in a usual sense, but rather as distributions in loop space. The momentum loops which 
we compute instead, are expected to be the usual numbers, with smooth momentum 
dependence beyond perturbation theory. 

Once again the large N matrix technology proved its power. Our matrix model is 
manifestly 0{d) invariant in any space dimension. Unlike the Eguchi-Kawai model 0, it 
has no matrix integrals, but rather algebraic equations. In a way, this is simpler, as we 
could use variational methods with some trial functions, dictated by physical intuition. 

The Fock space of our theory is remarkably small. In fact, this is the very similar to 
we would get in the Eguchi-Kawai model. This was the model oi d = 4 random matrices 
with some invariant probability measure. Treating this model by means of the operator 
expansion of the noncommutative probability theory we would get the similar Fock space. 
In four dimensions this a space of all sentences made from the four letter words - an 
infinite space, but still a discrete one. 

There is a significant difference though. The EK model was the lattice model, with 
inevitable large N phase transition. The operator expansion in this model would cor- 
respond to the strong coupling expansion of original lattice gauge theory. Here, on the 
contrary, the local limit is already taken. We are dealing with continuum theory, and we 
are expanding the Lorentz invariant amplitudes in the low energy-momentum limit. 

In terms of the Wilson loops, these amplitudes are given by path integrals over all 
paths, with Fourier measure (not the Wiener measure!). The absolute convergence of 
these integrals is to be achieved only at the nonperturbative level, due to the area law. 
At small momenta, the dominant paths would cover the area of the order of inverse string 
tension, which would lead to the momentum expansion coefficients proportional to inverse 
powers of the string tension. 

These are not the directly observable quantities, because we do not include the quark 
degrees of freedom at the endpoints of our string. Those are included later in the meson 
Hamiltonian we construct. There are two terms, the free particle one, plus the contribu- 
tion from the master field. This field is given by the gradient of the delta function of the 
X operator. The matrix elements of such master field between the trial states generated 
by the vertex operators are well defined and calculable. 

In the glueball sector we derived the linear wave equation which is supposed to describe 
the mass spectrum. Again, the coefficients of this equation are directly related to the 
coefficients of the operator expansion. Truncating this equation would result in the finite 
linear equation in the glueball sector. In a way, this is even simpler than the meson wave 
equation. 

We did not carry through these calculations, bat rather sketched the method. Working 
out the details and deriving these wave equations would be the most urgent next step. 
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An obvious question is how this agrees with the 't Hooft solution |T]] in two dimensions. 
I do not know the full answer to this question, but our meson Hamiltonian in the light 
cone gauge t = x+ could be compared to the 't Hooft's one. Clearly, the interaction term 
in the 't Hooft's Hamiltonian corresponds to the string momentum operator P_ which is 
related to our position field as discussed above. 

So, we could interpret the 't Hooft solution in our terms. The real question is how to 
get this solution from our equations. This would give us an insight for the solution of the 
4D problem. The momentum loop equation can be studied in the light cone gauge, where 
all the dynamics reduces to Coulomb "gluons" with propagators 4-- 

The most appropriate would be the mixed x+,p- representation, where the momentum 
loop equation must reproduce the Bethe-Salpeter equation for the ladder diagrams. As 
was pointed out by 't Hooft, the x+,p_ propagation goes along the strip with fixed width 
in p„ direction. It would be an interesting exercise to push forward this analogy. 

Regardless the 2D tests, the suggested approach could be tested in 4D by recovering 
the planar QCD graphs. The formal perturbation expansion follows from direct iterations 



of the momentum loop equation [ 12 1, assuming the bare vacuum as a zeroth approximation. 

Moreover, the perturbative QCD with running coupling must be recoverable in our 
representation at large p , or small X. This must be an interesting computation, making 
use of some WKB asymptotics. 

Finally, why not try to numerically compute the hadron spectrum, truncating the 
operator expansion and using variational approach like the one in[0]? This is a hard 



problem, but with the algebraic formulation we are getting closer to its solution. 
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A Functional Fourier Transform 

Let us parametrize the loops by the angle 6 which varies from to 27r. The functional 
measure for the path is defined according to the scalar product 

r rif) 
{A,B)=f-A{e)B{e) (140) 
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which diagonahzes in the discrete Fourier representation (not to be mixed with functional 

Fourier transform!) 

+00 

A{e) = Y.A^e'^'- A^^ = Al (141) 

—00 

+00 00 

(A, B) = Y, A^B_^ = AoBo + E a'X + a'X', < = V23iA,„ < = V2QA^ (142) 



The corresponding measure is given by an infinite product of the Euchdean measures 
for the imaginary and real parts of each Fourier component 

00 

z}g = rf'^gonA;An (143) 

1 

The orthogonality of Fourier transformation could now be explicitly checked, as 

f DCexp (t fdeCaie) (A^ie) - B^{e))] (i44) 

= J d'^Co n d%d''c: exp (2m (co (Ao - 5o) + E cL «, - b'J + c^ K ' K)\ ) 

= 5\A,-B,)X{5'K-K)5'K-K) 



Let us now check the parametric invariance 

6 ^ f{e)- /(27r) - /(O) = 27r; f\e) > (145) 

The functions x{6) and p{6) have zero dimension in a sense, that only their argument 
transforms 

x{e)^x{f{e));p{e)^p{f{e)) (i46) 

The functions x'{6) and p'{6) have dimension 1 

p'iO) ^ f'i9)p' if (9)) ; x'i9) -. f'i9)x' (/(^)) (147) 

The invariance of the measure is easy to check for infinitesimal reparametrization 

f{9) = 9 + e{9y, e(27r)=e(0) (148) 

which changes x and {x, x) as follows 

/d9 t d9 

e{9)2x^{9)x'^{9) = -i —e'{9)xl{9) (149) 

ZTX J Zn 

The corresponding Jacobian reduces to 

I- ld9e\9) = l (150) 
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in virtue of periodicity. 

The computation of the functional determinants can be performed in terms of discrete 
Fourier harmonics 

J = j Dx5\x{Q))e^^{^iF^, j dex^{e)x'M) = (151) 

oo „ 

n / rf'^<rf'x::exp {27,1 nP,, (x^)^ «),) . 

n=l-' 

This yields the infinite product 

00 -. 

J =\{{n'^ \deiF\y , (152) 

n=l 

which we define by means of the C, regularization 

J= lim^exp (-^ln(n'^ |detF|) [n'^ \<\eiF\)] (153) 

This yields 



\ n=l 



J = const |detF| '^^"^ = const ^J \([ei F\ (154) 

In odd dimensions this is zero, while in even dimensions it gives the pfafiian of F. 

The easiest way to obtain this result would be to note, that it comes from elimination 
of the zero mode. In the original form we could have locally rotated variables 

x^{e) ^ n^,x,i9), (155) 

to reduce the antisymmetric matrix to a canonical Jordan form 

- j Ff^^x^dx^ ^ ^ /^ y p^dqk (156) 

fc=i 

Then we could rescale the coordinates 

_^ Pk _^ Qk /-, p.„N 

Pk => -7^, Qk => -7f= (157) 

VJk VJk 

after which the matrix F would disappear from the exponential. Only the delta function 
would acquire the corresponding Jacobian 



6{x{0)) => H fkS{x{0)) = VdetF6{x{0)). (158) 



k=l 

This is the correct factor, so the remaining jacobian from rescaling of the variables in the 
measure should cancel. 

This is, indeed, so. Formally, the jacobian of this transformation is 

(detv^)"^^°\ (159) 
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where S(9) is the periodic delta function 

oo oo 

S{e) = ^ exp (z n^) = 1 + 2 ^ cos(n^). (160) 

— oo 1 

Its hniit at ^ — > can be defined by analytic regularization 

oo 

S^{9) = l + 2Y,n''' cosine). (161) 

1 

We find 

5„(0) = 1 + 2C(a) (162) 

which in the limit a ^ yields zero, as it should. This result means, that our measure 
with analytic regularization is scale invariant x{6) =^ Xx{6) in addition to parametric 
invariance x{9) ^ x{f{6)). 

B Numerical Approach 

For the numerical solution these relations could be globally fitted by some variational 
Anzatz for X^ with finite large A^, like it was done before [0 for the discrete version of 
the momentum loop equation. 

The discretization used in that work also corresponded to a polygon in momentum 
space, but the discrete equations were valid only up to the powerlike corrections at large 
number L of vertexes of this polygon. The infiuence of this discretization on the renor- 
mahzability of QCD is yet to be studied theoretically. It could be that any approximation 
of this kind to the singular integral equations destroys asymptotic freedom. 

The trial function was the Gaussian, 

M{ki,...kL)= f Dxexpl-- f xA W{x) exp (^-tY,kix{ti)) (163) 

\ uj=2ivn/T 

The discrete Fourier expansion at the loop was used. 

^ 1 

p{t)=Y.kie{t-ti) = - Y. e*"V-, P-c.=p::,Po = 0. (164) 

1=1 uj=2Tm/T 

The variational function Q{uj) was parametrized by eight parameters (see[|15|). 
The Regge slope in this model is given by 

ZTT 
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The discrepancy of the momentum loop equation was squared, analytically integrated 
over all momenta (with the Gaussian cutoff exp f — Y,t=i ^a^ ) ) ^^^ then minimized for 
L = 3, 4, .... The numerical limit L = oo was reached at L ~ 40. 

Numerically, the renormalization group behaviour 

was shown to agree with the variational solution of the discretized momentum loop equa- 
tion, but the accuracy in equation (about 0.1%) was still not sufficient to derive renor- 
malization group behaviour numerically (see[0 for more details). 



The important advantage of the present set of equations for the dual amplitudes is that 
no approximations were made. These are exact equations, therefore, the variational solu- 
tion with the general enough Anzatz must reproduce asymptotic freedom. For example, 
one may try the above WKB approximation with adjustable running coupling. 
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